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Abstract 

The problem of computing the one-dimensional configuration sums of the ABF model in 
regime III is mapped onto the problem of evaluating the grand-canonical partition function 
of a gas of charged particles obeying certain fermionic exclusion rules. We thus obtain a new 
fermionic method to compute the local height probabilities of the model. Combined with 
the original bosonic approach of Andrews, Baxter and Forrester, we obtain a new proof of 
(some of) Melzer's polynomial identities. In the infinite limit these identities yield Rogers- 
Ramanujan type identities for the Virasoro characters Xii (?) as conjectured by the Stony 
Brook group. As a result of our working the corner transfer matrix and thermodynamic Bethe 
Ansatz approaches to solvable lattice models are unified. 

Key words: Corner transfer matrices; Thermodynamic Bethe Ansatz; ABF model; One- 
dimensional Fermi gas; Roger s-Ramanuj an identities. 



1 Introduction 

The thermodynamic Bethe Ansatz (TBA) and the corner transfer matrix (CTM) method are two 
of the most fruitful techniques developed for studying solvable lattice models. The TBA approach, 
originating from the work of Yang and Yang Jl| and refined and extended in the work of, a.o., 
Takahashi 0, Takahashi and Suzuki [[|, and Bazhanov and Reshetikhin [|], [|, []], has been applied 
most successfully in studying the thermodynamic properties of the energy spectrum of the row-to- 
row transfer matrix. Also, in establishing connections between (off-)critical solvable lattice models 
and (perturbed) conformal field theories, the TBA method has proven to be extremely useful. The 
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CTM method in turn, was invented by Baxter (see ref. |7|] and references therein) and developed 
for computing order parameters of solvable models. 

At first sight, both methods seem rather distinct. Different types of transfer matrices are used 
for computing different quantities of interest. One key-step in both approaches is however similar .0 
Given a model system on a finite lattice of size m, one is naturally led to define a "truncated" or 
"finitized" system, which in the limit m —>■ oo becomes identical to the true model of interest. Let 
us describe this step in some detail as it turns out to be of unexpected importance. 

In using the Bethe Ansatz technique, one has to solve a coupled set of non-linear equations, 
known as the Bethe Ansatz equations (BAE). These BAE have many solutions corresponding to 
the different eigenvalues of the row-to-row transfer matrix. In solving the BAE numerically for small 
systemsizes, one observes that the roots of the BAE occur in so-called strings, a string being a set 
of complex numbers with equal real part but distinct imaginary parts. For large m, the imaginary 
positions of the roots within a string tend to a limit (fixed by the parameters in the model) and 
the only degrees of freedom left in the problem are the positions of the strings on the real axis, 
and the "occupation" number of a string of given type. For finite m the string-picture is violated 
by: i) The roots within a string do not exactly have equal real part, i.e., the strings are generally 
slightly curved, ii) The imaginary parts of the roots within a string deviate from their fixed m — > oo 
position. Hi) A number of roots does not fit the string picture at all, where it is to be remarked 
that the number of such roots divided by the total number of roots tends to zero in the large m 
limit .0 

Based on the above observations, one usually proceeds by making a string hypothesis, which 
amounts to assuming that all solutions to the BAE consist of sets of strings, each string being taken 
from a set of allowed string types. Of course the string hypothesis, if correct at all, will only describe 
the spectrum of the row-to-row transfer matrix in the m — > oo limit, and the set of "reduced" Bethe 
equations obtained after substitution of the strings makes sense in the thermodynamic limit only. 
It is nevertheless tempting to view the reduced equations as those of some ideal solvable model for 
which the string hypothesis holds true even for finite m. We call this virtual model the finitized 
model. Of course, the finitized model should not be confused with the original model on a finite 
lattice, but should be viewed as a model on a finite lattice inheriting all simplifying features of the 
original model when the lattice size is taken to infinity. 

A similar finitization occurs in the CTM calculations. Provided the corner transfer matrices 
commute, Baxter showed that all eigenvalues of the (diagonalized) CTM's of a solvable model 
take a special exponential form, which can be computed explicitly by considering the ordered low- 
temperature limit. Since commutativity of the CTM's holds only in the infinite lattice limit, we can 
only diagonalize the CTM's in this limit. Let us think of the exponential form of the eigenvalues as 
being at the same level as the string hypothesis. Again we can substitute this form in the relevant 
finite m equations and think of the equations thus obtained as those of an ideal solvable model for 
which the CTM's commute for finite systemsize. We again call this the finitized model. 

The aim of this paper is to show that given a particular solvable model the finitized models 
obtained from the TBA and from the CTM approach are the same. To be precise, we show that the 

1 Throughout this paper we only refer to models that satisfy the Yang-Baxter equation with difference property 
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2 Recently solutions to the Bethe equations for the 3-state Potts model have been found which are neither of string 
type nor thermodynamically irrelevant. For the full details on such non-string solutions see ref. fa. 



finite m TBA and finite m CTM calculations for the finitized model are mathematically equivalent. 
Letting m to infinity, this establishes the claim that the true (infinite m) TBA and CTM methods 
are equivalent as well. We proof the above assertion for what is probably the most generic example 
of a solvable lattice model, the (r — l)-state Andrews-Baxter- Forrester (ABF) model 0. To do 
so we present a new fermionic method for computing the one-dimensional configuration sums of 
the ABF model. This method is based on a mapping of the configuration sums onto the grand- 
canonical partition function of a one-dimensional gas of charged fermions. Up to terminology, this 
gas is equivalent to the system of strings and holes describing the TBA. 

The remaining sections of this paper are organized as follows. In section |2| we briefly review 
the definition of the ABF model and sketch the bosonic approach for computing the local height 
probabilities as followed by Andrews et al. Then, in section || we present the main ideas of our 
fermionic method by computing the simplest possible configuration sum, using the Fermi gas. In 



the subsequent section we discuss Melzer's polynomial identities and their limiting Rogers 



Ramanujan type identities ||11| , which are proven by combining the bosonic and fermionic methods. 
In section f| we also discuss the relevant differences between our approach to proving Melzer's 
identities and the original proof given by Berkovich [T2|] . In section |^, we discuss the equivalence 
between CTM and TBA as follows from the working of section |3], and we end with a summary and 
discussion in section |5|. 

In this paper we have limited ourselves to the fermionic evaluation of just a single configuration 
sum. The computation of all configuration sums will be the content of a subsequent paper, referred 
to as part II. In this part II, we show that two additional boundary particles with fractional charges 
have to be introduced on top of the Fermi gas to treat the most general configuration sums. A brief 



explanation of the origin of these boundary particles is given in section 3.3 



Before we start reviewing the aspects of the ABF model relevant to this paper, several final 
introductory remarks have to be made to put some of the ideas and results of this paper in proper 
context. First of all, we have to mention the remarkable series of papers written by the Stony Brook 
group [11,13-15]. Based on a study of the completeness of states of the Bethe Ansatz solution for 



the 3-state Potts model [13, M, large classes of new expressions for conformal characters where 



conjectured Hl6 , 11 1. Because of exclusion rules labelling the Bethe states, these character expres- 
sions were called fermionic. Combining the fermionic forms with the well-known Rocha-Caridi type 
bosonic character expressions, led to conjectures of many new identities similar to those of Rogers 
and Ramanujan. Second, In ref. [17]| Jimbo et al. pointed out that computing one-dimensional con- 
figuration sums using the recursion method of ABF leads to bosonic character expressions. This, 
together with the fermionic forms originating from the Bethe Ansatz approach indicates a deep 
connection between (T)BA and CTM methods. Signs of an even stronger connection were found 
by Melzer [ID|, who conjectured that finitizations of the fermionic character expressions, originating 
from TBA for the finitized ABF model, equal finitizations of the bosonic character expressions orig- 
inating from CTM for the finitized ABF model. It has especially been this last result and its many 
generalizations to other solvable models [18-26], that has motivated our attempt to unify TBA and 
CTM approaches to solvable models. 



2 The ABF model 



This section serves as a brief introduction to the ABF model, and, apart from the model definition, 
reviews part of the calculation of the order parameters for regime III as carried out in ref. ||. 



2.1 Definition of the model 



The ABF model is a restricted solid-on-solid (RSOS) model defined on the square lattice C. Each 
site % of C carries a spin or height variable crj, which can take any of the values 1, 2, . . . , r — 1. A 
restriction is imposed on each pair of neighbouring sites % and j by |<7, — <Jj |=1. The Boltzmann 
weight of a configuration of heights on the lattice is computed as the product over elementary 
face-weights, with nonzero faces parametrized as 
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Here we have used the elliptic theta function with argument u and arbitrary but fixed nome p, 
\P\ < 1, 

oo 

tfi (u) =2j9 1/8 sin M n (l -2p n cos2M + p 2n ) (l-p n ). (2.2) 

n=l 

The variable w in (|2.1| ) is the usual spectral parameter, and A is the crossing parameter, fixed by 
the number of states of the model as 

A = -. (2.3) 
r 

In ref. the weights W were shown to solve the Yang-Baxter equation 0, and hence quantities 
like the free energy and order parameters can be computed exactly. In their calculations, ABF 
distinguished four different physical regimes, but here we restrict ourselves to regime III, given by 

regime III: 0<w<A, < p < 1. (2.4) 



2.2 Computation of the local height probabilities 

As a first step in calculating the order parameters of the ABF model (see ref. |]27| for the actual 
definitions), one has to compute the so-called local height probabilities P bc (a). Here P bc (a) is the 
probability that a spin of the (infinite) lattice takes the value a provided that the model is in 
a phase indexed by b and c. In their computation of P bc (a) ABF start with a finite lattice £ 
with size measured by m, and consider the probability P^(a) that the center spin of C takes the 
value a, provided the boundary spins on C are fixed according to a groundstate labelled by b, c. 
For regime III there are 2r — 4 antiferromagnetic groundstates, corresponding to all spins on the 
even sublattice taking the value b and all spins on the odd sublattice taking the values c = b ± 1, 
b = (3 =f l)/2, (5 =f l)/2, . . . , (2r — 3 =F l)/2. We can thus write 



P bc (a) = Jim P%(a), (2.5) 



with b, c taking any of the above groundstate values. At the same time ABF showed, using the 
corner transfer matrix method, that 



pbc( n \ _ H„. E(x a 1 x r )X m (a,b, c; x 2 ) 

with the one- dimensional configuration sums defined as 

X m (a,6,c;g) = £ gE^i*!^-^!/* aQ = a> ^ = 6j am+1 = c . ( 2 .7) 

0"l,...,(T m _l = l 

The nonstandard labelling of the ex's from to m + 1 is chosen for later convenience. The variable 
x in (|2.6| ) relates to the nome p of the Boltzmann weights as 

x = e -**/(re) ) with p = e~ e . (2.8) 
The function E in (|2.6|) is the elliptic function 

oo 

E(z,p) = l[(l-p n - 1 z)(l-p n z- l )(l-p n ), (2.9) 

n=l 

for all z,p E <£, \p\ < 1. 

We remark here that in deriving ( |2.6| ), one has to assume commuting CTM's, and hence that the 
size of the lattice is infinite. Therefore, the parameter m in ( |2.6| ) is not the same as the systemsize 
m in Q2.5| ). As described in the introduction, we now define our finitized model by identifying the 
m's in (|2.6|) and (p.5|). In other words, we define the finitized model through P^{p) as 
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E a =i E(x a , x r )X m (a, b, c; x 1 ) 



The problem now is to perform the sum in ( |2.7| ) to obtain manageable closed form expressions 
for the configuration sums X m (a, b, c; q). In the remainder of this paper we will give two approaches 
to the solutions of this problem. The first is the method originally employed by Andrews, Baxter 
and Forrester . It results in expressions for the configuration sums as the difference of two finite 
g-series, and hence, adopting the terminology introduced by the Stony Brook group [0, PH , we 
call this method bosonic. The second method is new, and amounts to interpreting the sum in ( |2.7|) 
as the grand-canonical partition function of a one-dimensional gas of charged particles, obeying 
certain exclusion rules. This time the resulting expressions involve just a single finite g-series. 
Again following ref. [16], |H , we call the second method fermionic|] 



3 We remind the reader that Baxter's original method of computing P bc (a) for the hard hexagon model also led 
to infinite fermionic g-series |E8[ . Our method can thus be seen as a generalization of the working of ref. p8| . 



2.3 Bosonic evaluation of the configuration sums 

Let us begin with a brief reminder of the original bosonic approach of Andrews, Baxter and Forrester, 
as detailed in ref. M. 

First, notice that the X m (a, b, c; q) can be (re)defined in terms of the following recurrence rela- 
tions: 

X m (a, b, b ± 1) = X m „x(a, b ± q, b) + q m/2 X m ^(a, 6^1,1), 3 < 2b ± 1 < 2r - 3, (2.11) 
subject to the initial and boundary conditions 

X (a,b,c) = 8 a> b8 Ctb ±x, 



X m (a, 0, 1) = X m (a, r, r - 1) = 0. 
Then, introducing the Gaussian or g-binomials as 
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with (q) m = Ilfcli(l ~ Q k ) for m > and (q)o = 1 , we have the following theorem ||: 

Theorem: For m > 0, 1 < a, b, c < r — 1, c = 6 ± 1, m + a — b E 2> , the one-dimensional 
configuration sums read 



X m (a,b,c) = g(-<»(-)/ 4 £ <7 
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For the proof of this result we refer the reader to ref. [[J. Let us however remark that it follows in 
straightforward manner substituting ( |2.14 ) into ( [2.11 ) and using the elementary recurrences [2S 
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Also the conditions ( |2.12 ) are readily checked for the expression ( |2.14| ). 

Before concluding this section we note that the X m (a, b, c;q) can be viewed as finitizations of 
the Rocha-Caridi expressions for the (normalized) Virasoro characters X(6+c-i)/2a(?) 0> 

lim q-^ a -^ 4 X m (a,b,c;q) 
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3 Fermionic evaluation of the configuration sums 



We now come to the main part of this paper, the alternative method to evaluate the sum Q2.7Q . 
As mentioned earlier, we can interpret this sum as the grand-canonical partition function of a one- 
dimensional gas of charged particles. Once this identification has been made, the problem of course 
remains to actually compute this partition function. For pedagogical reasons this is in fact what we 
do first, that is, we start by defining our system of particles and compute its generating function. 
Only afterwards we then perform the actual mapping from the one- dimensional configuration sums 
onto the Fermi gas partition function. To keep things as simple as possible we restrict ourselves 
to the evaluation of X m (l,l,2). The case of arbitrary X m (a, b, c), which involves the additional 
concept of boundary particles with fractional charge, will be treated in part II of this paper. 

3.1 A one-dimensional Fermi gas 

We consider a system of charged particles occupying a one-dimensional lattice of m + 1 sites, with 
sites labelled from to m, m being even. The number of particles with charge j is nj (j = 1, . . . , r— 2) 
and a particle of charge j has diameter 2j. The lattice is completely occupied with particles, yielding 
the completeness relation 

r-2 

A typical configuration of particles is shown in figure [l|, where we have drawn a particle of charge 
j as a triangle with height j. 

To describe the allowed motion of the particles, we have the following set of rules: 

Rl Hard-core repulsion between particles of equal charge. Therefore, given an "initial" configu- 
ration, the order of the particles with charge j remains fixed. For obvious reasons we can 
therefore think of the particles as being fermions. 

R2 Two particles with charge j and k, j ^ k can penetrate each other, and eventually, exchange 
position. An example of the motion of a particle with charge k through a particle with charge 
j > k is shown in figure |2|. We note that the diameter of the charge complex in this figure 
remains 2(j + k) = 2j + 2k throughout the motion and hence that we obey charge conservation. 
We also note that both the position and the number of peaks of the complex remain fixed. 

R3 We identify configurations that lead to the same profile as being the same. Hence the charge 
configurations c and c' in figure || are identical. (Although Rl is a special case of R3, we wish 
to state it as a separate rule.) 

With the above rules it is clear that a configuration as drawn in figure [3] is not generic, since all 
particles are properly separated and no charge complexes occur. A more general situation, which 
can be obtained from the configuration of figure [I] by moving around the particles following the 
rules of motion R1-R3, is shown in figure |3|. 

We now come to a crucial observation about our Fermi gas. Fix the number of particles of 
given charge, i.e., fix n 1; n 2 , . . . , n r _ 2 , and call the configuration with all particles being separated, 
and all particles of charge j being positioned to the left of all particles of charge k < j, a minimal 



configuration.^ Then all other configurations with the same particle content, i.e., with the same 
ni, n 2 , . . . , n r „ 2 can be obtained by the following motion starting from the minimal configuration: 

Ml Keep the particles with charge r — 2 fixed. 

M2 Move the left-most particle of charge r — 3, denoted p r -3,i to the left through some (and possibly 
all) of the particles of charge r — 2. It can either end up as a charge complex or as a separate 
particle. 

M3 Move the second-left-most particle of charge r — 3 (denoted p r - 3,2) to the left. Of course, this 
time the movement to the left is bounded by the position of particle p r -3,i due to the fermionic 
exclusion rule Rl. 

M4 Continue to move all particles p r -3, n , n = 1,2,..., n r _ 3 to the left, one at the time and p r _ 3in+ i 
after p r -3, n - Each time the motion of particle p r _ 3jn+ i is bounded to the left by the position 
of particle p r -3, n - 

M5 After having moved the particles of charge r — 3, start a similar procedure for the particles of 
charge r — 4. The left-most particle p r -4,i first and right-most particle p r -4,n r _ 4 last. 

M6 Continue this process till all particles have moved to the left, always obeying the order: p^ n 
before p^ n i for n < n' and Pj >n before p^ n i for j > k. Each time the particle Pj >n +i is bounded 
to the left by pj, n , with pj \ bounded by the origin of the lattice. 

In figure ^ we have shown how to obtain the configuration of figure |3| by consecutively moving the 
particles of the corresponding minimal configuration leftwards. To proof that we indeed can obtain 
all allowed configurations by starting from a minimal configuration carrying out the above steps, 
we first note that given a configuration, e.g., that of figure |3|, we can uniquely determine its particle 
content in the following way: 

• First locate all charge complexes, regarding a separate particle as a complex (though trivial) 
as well. This amounts to marking all coordinates on the lattice that have zero height. In the 
case of the configuration of figure ||, the first complex ranges from to 6, the second from 6 
to 22 and the third from 22 to 34 = m. 

• For each separate complex we now determine its particle content as follows: 

— Draw a dashed line along the lattice from the origin (left-most point) of the complex to 
its endpoint (right-most point). We call this the zeroth baseline of the complex. 

The result of this trivial step is shown in figure ^a for a typical charge complex. 

— Start from the highest point of the complex. If there is more than one highest point, 
start from the left-most highest point. Move down to the right of the peak along the 
contour of the complex till its endpoint. (The point were the next complex starts.) If 
a local minimum is reached, i.e., the contour of the complex starts going up again, we 
draw a dashed line from this local minimum to the right, until we cross the contour of the 

4 This terminology, suggested to the author by O. Foda, stems from a partition theoretic interpretation of minimal 
configurations. 



complex. At that point we move further down along the contour. If another minimum 
occurs we repeat the above, et cetera. 

We do exactly the same, now starting from the same left-most highest point but moving 
down to the left. If a local minimum is reached we draw a dashed line to the left and 
continue our movement down when the dashed line intersects the contour of the complex. 

The result of the above procedure is shown in figure |5|b. The total of dashed line segments 
drawn in this step is called the first baseline. 

— We now view the first baseline as the zeroth baseline of smaller complexes. Let there be 
L such smaller complexes. We then divide the first baseline into L pieces by marking 
the begin- and endpoint of these complexes by arrows pointing to the left and right, 
respectively, see figure ||b. 

Either these complexes consist of a single particle (complexes with one maximum and 
consequently no local minima), or they are again complexes in the true sense of the word. 
In the latter case, we repeat the procedure of drawing dashed lines, again starting from 
the left-most highest point. 

The result of this third round of drawing dashed lines is shown in figure ^c. The total 
of dashed line drawn in this step is called the second baseline. Again this baseline is 
divided into pieces, marking yet smaller complexes. 

— We repeat the above procedure of drawing and cutting up higher order baselines itera- 
tively. That is, we view the j-th baseline as the zeroth baseline of yet smaller objects, 
we draw arrows accordingly and, if some of these smaller complexes have more than one 
maximum, we start drawing the (j + l)-th baseline. 

Since the highest point of the complex is always less than r — 1, the iteration terminates 
within r — 2 steps (including the drawing of the zeroth baseline). In the case of figure |5], 
we are done in four iterative steps, figure [|d being the final result. 

— Given the complex with all the baselines we can read of its particle content straight 
away. First of all, according to R2, each peak corresponds to a particle. To determine its 
charge we move down vertically from its peak till me meet the zeroth baseline relative to 
this particle, i.e., we move down till we first intersect a baseline. The height of the peak 
minus the height of the baseline at the intersection point is its charge. In the example 
of figure ||d, we thus get (ni, n,2, ■ ■ .) = (2, 5, 1, 0, 1, 0, 0, . . .). 

• Clearly we now determine the content of a given configuration as the sum of the particle 
content of all its complexes. 

We note that in solving the problem of determining the particle content of a configuration we have 
actually localized each individual particle. Hence we can start moving the particles to the right to 
obtain the minimal configuration. We perform this motion in exactly the opposite order as described 
under M1-M6. Since our rules of motion are reversible (left-right invariant) we have consequently 
proven our assertion that each configuration can be obtained from a minimal configuration following 
the previously described rules. 

A final ingredient needed in the Fermi gas is the actual Boltzmann weight of a given charge 
configuration. Given this we can try to compute the partition function. Before we describe this 
missing piece of information, let us first compute the total number of configurations keeping the 



particle content fixed, as this naturally leads us to introduce some final notation needed. Consider 
the particle of charge j as shown in figure |6|. We denote the left-most point of the particle its origin, 
and its right-most point its endpoint. All points 1, 2, . . . , 2j — 1 will be referred to as interior points 
of the particle. Let us now ask the question in how many ways k) this particle and the particle 
with charge k < j shown in the same figure, can form a charge complex. Clearly, (see also figure 0), 
we can place the origin of the particle with charge k at the following interior points of the particle 
of charge j: 

2j-l,2j,...,j + k. (3.2) 

Similarly, we can place the endpoint of the particle at the following interior points: 

1,2,..., j-k. (3.3) 

If we now remark that, according to rule R3, placing the origin at j + k gives the same configuration 
as placing the endpoint at j — k, (see figure 0c and c') we have a total of 

I(j,k) = 2j-2k-l (3.4) 

ways to form a complex of charge j + k. 

This prepares us to answer the combinatorial question raised above. The number of ways to 
move the particles with charge r — 3 to the left, starting from the minimal configuration, is 



'ff f \ . . ^ff 1 = ( 2nr - 2 + nr - 3 ] , (3.5) 



fci=0 k 2 =0 k„ r ,=0 



n r _ 3 



since I(r — 3,r — 4) = 1, and since we can position a particle not only in the interior of a larger 
particle as a charge complex, but also in between two larger particles as a separate unbound particle. 
Similarly, the number of ways to move the particles with charge r — 4 to the left is 

E E • • • E 1 = > (3-6) 

fc 1= k 2 =0 k nr _ 4 =0 \ '^-4 J 

since J(r — 2, r — 4) = 3 and J(r — 3, r — 4) = 1. In the general case with particles of charge j we 
obtain 



2(r-2-j)n r _ 2 +...+4n i+ 2+2?ij + i fc 1 



1=0 k n .=0 

2(r — 2 - j)n r _ 2 + . . . + 4n,- +2 + 2nj + i + ' 



fei=0 fc 2 =0 k n .=0 
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(3.7) 



as J(r — 2, j) = 2(r — 2 — j) — 1, ... , /(j + 1, j) = 1. Of course the above calculation is rather clumsy, 
and we could have written the binomial answers straight away. As we later need the g-analogues of 
(|3.5| )-(3~7) we nevertheless thought it instructive to present the result in the above written form. 

Collecting the above results we get for the total number of configurations Z with fixed particle 
content n\,n%. . . n r _ 2 , 

Z(n u u 2 ,..., n r - 2 ) = if ('tt* - 2 - i - *K-« + n,\ (J _ g) 

i=i V n i / 



In turn summing over different particle content, computing the number of configurations S m in a 
grand-canonical setting, we get 



Z(m,n 2 , ■ ■ ■ ,n r _ 2 ), (3.9) 

ni,n2,...,n r _2>0 



where the prime over the sum indicates the restriction ( |3.1|) , i.e., we sum over particle numbers 
keeping the size of the system fixed. 

This last result suggests in fact to eliminate one degree of freedom. To do so we introduce the 
variables rrij as the number of antiparticles of charge —j, 

m i = 2n i+ i + 4n i+2 + ... + 2(r-2-j)n r _ 2 , j = 1, . . . , r - 2. (3.10) 

From this we see that antiparticles of charge — (r — 2) are absent, and hence we can eliminate the 
dependence on the particles of charge r — 2. Using the vector notation m = (m 1 ,m 2 , . . . ,m r _ 3 ) T , 
n = (n 1; n 2 , . . . , n r _ 3 ) T and = and denoting the incidence matrix of the A r _ 3 Dynkin 
diagram as X, we find 

m + n = - (X fh + m ei) , (3-H) 

and 



Here the sum is over all (integer) solutions to the constraint equation ( |3.11|) , where the prime 



denotes the restriction to even antiparticle numbers rrij, as follows from ( [3.10|) . Equation ( [3. 11] ), 
relating the variables rij to the new variables rrij, originates from the work on fermionic sums by 
Berkovich [12]. As we will see in section |5|, it also appears in the context of TBA. 



3.2 Computation of the partition function 

We now finally come to the actual definition of the Boltzmann weight of a charge configuration, 
and to the calculation of the grand-canonical partition function. 

Let us return to the charge configuration of figure ^. Clearly, at each (integer) point x (x = 
1, ... ,m — 1) of the lattice there are four possibilities for the contour of the configuration.^ We 
either have a straight line going up or down, or there is a cusp corresponding to a (local) minimum 
or maximum. We now assign an energy E(x) to each site x of the lattice by 

n _ J ^xlogq if the contour at position x is a straight line , , 

I if the contour at position x is a cusp. 

Clearly, with this definition the groundstate of the model is given by the state with only particles 
of charge 1. For given, fixed particle content, the energy is minimized by ordering the particles to 
form a minimal configuration. 



To be precise we should of course say "contour in a small neighbourhood of position x." 



The remainder of this subsection is devoted to the computation of the grand-canonical partition 
function of our Fermi gas, 



= E Z(ni,...,nr- 2 ;q), (3.14) 

particle content 

with the partition function of fixed particle content given by 

m— 1 

-P E E{k) 

Z(n 1 ,...,n r _ 2 ;q)= e k=1 ( 3 - 15 ) 

charge configurations 

To compute Z we follow a procedure similar to that employed in the section f>.l\ First we consider 
the minimal configuration, and from that we obtain all other configurations with equal content by 
carrying out the steps M1-M6. The only extra input will be that changing the position of a particle 
changes the energy of a configuration. 

To compute the energy of a minimal configuration we proceed as follows. The energy of a particle 
with charge j, with its origin at position x is 

1 2j-l 

-f3E(j,x) = -logq £ (k + x) = (j + x)(j-l)\ogq (3.16) 

1 k=i, k^j 

Hence, for the minimal configuration with content n 1; . . . , n r „ 2 we compute 

r-2 «j / r-2 

-f3E min = \ogqJ2J2 E [j,2j(e-l) + 2 Yl kn k 

3=1 1=1 \ k=j+l 
r-2 / r-2 \ 

= loggE0' -1 ) n iU n j + 2 E kn k) ( 3 - 17 ) 

3=1 \ k=j+l J 

r— 3 i — 3 

= logg^E^'+i^^+i' 

3=1 k=l 

with matrix A defined as 

i>i (3 - 18) 

Writing this in favour of the occupation numbers of antiparticles rrij defined in ( |3.10|) , we can 
simplify to 

e -PE min = q \m T Cm^ (3ig) 

where C denotes the Cartan matrix of the simply- laced Lie algebra A r _ 3 . 

As before, all other states with the same content n\, . . . ,n r _2 are obtained from the minimal 
configuration by moving particles to the left following M1-M6. The combinatorics of this rearranging 
of particles has been considered in section |3.1| . What has not yet been worked out is the increase 
of energy associated with the leftward motion of a particle. Since we move one particle at a time, 
we only have to calculate the energy increase of the movement of a particle of charge k, relative to 



a charge complex with charge j > k positioned immediately to the left of this particle. Possibly 
the simplest way to go about is to first view the complex and the particle as one single complex 
with charge j + k. Suppose that relative to the origin of this complex we have (local) minima and 
maxima of the contour at < x\ < y\ < X2 < JJ2 < ■ ■ ■ < %n-i < Un-i < x n < 2k + 2j, labelling 
minima by y and maxima by x. Note that the simplest possible scenario; the complex of charge j 
consisting of a single particle, corresponds to n = 2. The total weight of the complex relative to its 
own origin is thus 



2j+2fc-l 

E '/ 

t = \ 

- 7^ -^ii V : 



k/2 



(3.20) 



Now start moving the particle with charge k one step to the left, like shown in figure [|. This results 
in a decrease of both y n -i and x n by 1. The weight increase of this step is hence a factor 



q 



(y n -!+x n +2)/2 



(j/„_i+x„)/2 



Q 



q. 



(3.21) 



Taking another step to the left, again leads to a decrease of y n -i and x n by 1 and we pick up another 
factor q. This process can be repeated until y n -i — % n -i — x n — y n -i, corresponding to the situation 
that the maxima at x n and x n -\ have equal height (see figure ^|c). From then on moving the particle 
a single step leads to a decrease of x n _i and y n -\ by 1. However, this still leads to an increase by 
a single factor q, and we are led to conclude that each step of particle k trough the complex yields 
an extra factor q. 

With this result we can simply carry out the g-analogue of the previous subsection. To do so we 
repeatedly make use of the identity pi] 



N 



fc 1= 0fe 2 =0 k M =0 



N + M 
M 



(3.22) 



First we start moving the particles of charge r — 3 by ( |3.22|) this gives a Boltzmann factor 

^'n r _3 — 1 



2n r -2 fei 

E E 

fc 1= fc 2 =0 



n r -3 



n r -3 + m r _ 3 
m r __ 3 



(3.23) 



Similarly, moving the particles with charge r — 4 to the left results in 

4n r _2+2n r _3 ki k„ r _ 



E • • • E q kl+k2 + --- + kn r~4 

fe 1= fe 2 =0 k„_,=Q 



4n r _2 + 2n r _ 3 + n r _ 4 
ra r _4 



n r _4 + m r _4 
m r -4 



.(3.24) 



J q 



In the general case of particles with charge j we end up with the contribution 



k nj -l 



2(r-2-j)n r _2+...+4n j+ 2+2nj + i ^ 

EE g fcl+fc2+ '" +fc " 
fci=0 fc 2 =0 



2(r 



k nj =0 



j)n r _ 2 + ■■■+ 4n j+2 + 2n j+ i + nj 
n 1 





rij + nij 


- q 


. m 5 . 



(3.25) 



From the results ( |3.23| )- (|3^25 ) and ( p. 19 ) we obtain an expression for the partition function Z 

as 



Z(n 1 ,n 2 ,...,n r _ 2 ;q) = q~ 



rh C rh 



r-3 

n 



rrij + rij 



(3.26) 



with equations ( |3.10| ) and ( |3.11| ) relating the antiparticle numbers rrij to the particle numbers rij. 
We can finally end this section by stating the result for the grand-canonical partition function, 



particle content 



rh C rh 



r-3 

n 



rrij 



(3.27) 



J q 



where, as in (|3.12|), the second sum is over all solutions to (|3.11|) with each entry of rh being even. 



3.3 Mapping of configuration sums onto Fermi gas partition function 



To show that the result (|3.27| ) is indeed the evaluation of ( |2.7| ) with a = b = 1 and c = 2, we fix the 



spins Oj G {1, 2, . . . , r — 1} such that \<Jj+i — <Jj\ = l for all j, and such that gq — a m = 1, cr m +i = 2. 
We call such a sequence (o"o, . . . , cr m+1 ) admissible. Next we plot Oj as a function of j and 
interpolate between (j,<Jj) and (j + l,a\, +1 ) by a straight line. We call the graph thus obtained 
the contour of (o"o, . . . , u m +i). A typical contour of an admissible sequence of sigma's is shown in 
figure 0. From the definition ( |2.7|) of the configuration sums we see that a straight line segment 
through (j,(Tj), corresponding to \(Tj-i — = 2, yields a factor q^ 2 . The total weight of an 

admissible sequence can therefore be written as 

m— 1 

e i =1 (3.28) 



where the energy function E(J) is precisely that of the Fermi gas given by equation ( |3.13| ). Since 



summing over all admissible sequences corresponds to summing over all contours, which in turn 
corresponds to summing over all charge configuration, we have established the desired equivalence. 

We note that particle-like interpretations of admissible contours (or paths) have been previously 
formulated for the minimal conformal series M(2, k) in refs. [31-33]. 

The above equivalence also explains the defining rules Rl and R3. From ( |2.7| ) we have to count 
each admissible sequence of spins just once. However, interchanging two identical particles leads to 
the same contour, and would therefore lead to multiple counting of one and the same admissible 
sequence. Similarly we have to impose R3 to avoid overcounting of admissible sequences. 

At the same time the mapping explains why the particular choice of X m (l, 1,2) is simplest to 
evaluate using the fermionic technique. If for e.g., we take a > 1, we have to introduce a particle, 
which, in separated unbound form, would have a profile of a single straight line going down. This 
boundary particle has a (fractional) charge of (a — l)/2. Similarly, when b = c — 1 > 1 we need 
a boundary particle of charge b/2 with profile of a straight line going up. These particles are in 
contrast with the particles introduced so far, which have origin and endpoint of equal height. Since 
the motion of a (bulk)particle through a boundary particle is very different from that defined in R2 
on page 0, the general X m (a,b,c) case becomes quite complicated and technical. For this reason 
the computation of all one- dimensional configuration sums will be deferred to part II. 



4 Polynomial and Rogers— Ramanuj an type identities 



Since the bosonic and the fermionic method to evaluate the configuration sum X m (l,l,2) yield 
rather different type of expressions, we can of course combine results to find some nontrivial identi- 
ties. Setting a = b = 1 and c = 2 in (|2.14|) , substituting fl3.11|) into (|3.27| ), and replacing m by 2m, 
yields 

j(Im)j + mSj i 



E 



l?n C ' fn 



mi , 



me22' 



-3>0 
-3 



r-3 

n 

3=1 L 



E 



r(r-l)j 2 +j 



2m 



m 



l)j 2 +(2r-l)j+l 



??? 



2m 
- 1 - 



(4.1) 



These identities are precisely the polynomial identities for X 2m (l, 1, 2) as conjectured by Melzer in 
ref. | r0|| . A different proof of ( ^4.1| ) has been published by Berkovich [Q. In his proof Berkovich 
has shown that a larger class of fermionic expressions, corresponding to all X m (l,b, b ± 1) solves 



the same recurrences ( 2.11|) as the bosonic expressions (|2.14 ) of ABF. To establish this, elegant 
generalizations of the decompositions (|2.15|) for the g-binomials were used. One of the great merits 
of the recursive proof of ref. |12[ is that it proves a whole set of fermionic expressions at the same 
time, and that complications with the boundary for general X m (a, b, c) are avoided | 3~4]| . On the other 
hand, the proof as given in this paper has the advantage of providing detailed information relating 
BA solutions (see next section), admissible sequences of spins occurring in the CTM calculation, 
charge configurations of the Fermi gas and actual terms within the sum on the left-hand side of 



l/(q)h and lim r , 



2m 



In taking the m — > oo limit of ( |4.1| ), using lim 
the following Rogers-Ramanujan type identities for the X[ r ^ 1,r \q) Virasoro characters arise: 

2 rh T C fn r-3 

e q -r, — n 



V(ff)c 



mi , 



,,m r _3>m r _2E 
mG2Z r - 3 



mi 



rrii 



.1 + m j+lj 



rrij 



1 oo 
1 E 



r(r-l)j 2 +(2r- 



(4.2) 



These identities were first conjectured by the Stony Brook group in ref. \.\ I 



To obtain Rogers-Ramanujan type identities for the branching functions of the Z r _2 parafermion 
conformal field theories of ref. |55| , we replace q by 1/q in ( [4.1| ) and again let m — ► oo. This results 
in 

,m T C _1 fn 

E 

mi,...,m r _3>0 
{C- x rh) r - Z e7Z. 

2 



(?)mi •••(?) 



m r _ 3 



(r-3)(r-2) 
q 24r 

[(^)oo] r 3 adQweW 



E E s s n( 



w) q 



a — 



rp — (r — l)w(p) 



r(r — 1) 



(4.3) 



with Q the root lattice, W the Weyl group and p the Weyl vector of A r _ 3 . We note that the left-hand 
side of these identities are the expressions for the branching functions obtained by Lepowsky and 
Prime To derive the above we have used the level-rank equivalence between the (r — l)-state 
ABF model in regime III and the level-2 A r _3 Jimbo-Miwa-Okado model in regime II |l7fl , to rewrite 
the SU(2) type right-hand side of (|4.1|) into an SU(r — 2) type form |^3J (see also the discussion in 
ref. f25|). For r = 4 we note that (4.2) and (|4.3| ) coincide. 



5 CTM versus TBA 

The reader familiar with the TBA computations for the ABF model of Bazhanov and Reshetikhin 
(BR) [[| , will have noticed the many similarities between the TBA and the Fermi gas. This section 
serves to indeed show the mathematical equivalence between the fermionic CTM calculations of 
section [3] and the TBA calculations of ref. [|J . 

First let us recall the Bethe Ansatz equations for the ABF model @, || 



UJ 



Here to = exp(i£A), £ = l,...,r — 1, A = n/r and each set of roots {cti, . . . , a m / 2 } yields an 
eigenvalue A(-u) of the row-to-row transfer matrix. 

Based on exact information for the cases r = 3 (trivial 2-state model) and r = 4 (Ising model), 
and on a numerical investigation for other r values, BR formulated the following string hypothesis: 
For m — > oo all solutions to the BAE consist of sets of strings, with the allowed string types given 
by 



U) = a U)_ll( J + 1 _2^ fi=l,...,j; j = l,...,r-2, (5.2) 



a 



" 2r 



where /i labels the j roots within a string and j the different types of strings. The real variable a^' 
is the string center. 

As described in the introduction, we define our finitized model by assuming the string hypothesis 
to be correct for finite m. Let nj denote the number of strings of type j. Then, since the total 
number of roots in a solution equals |m, we have the completeness relation 

r-2 

Y.3 n i = \ m - ( 5 - 3 ) 

We note this is exactly the completeness relation ( |3.1| ) for the particle numbers of the Fermi gas. 
Substituting the hypothesis into the finite m BAE, a set of reduced equations is obtained which 



only involve the string centers. Carrying out the standard TBA procedures 0, [37], the reduced 
equations yield the following set of constraint equations for the number of strings rij and the number 
of holes rrij (counting the missing Takahashi numbers): 

r-2 

m bj = m j + Bj,k n k , (5.4) 

k=l 



with 



2(r-j)k 



3>k 
j <k 



(5.5) 



Following BR we set j = r — 2 in to find m r _2 = 0. Now eliminating n r _2 yields precisely the 
constraint equation (|3.11|) , and we conclude that we can identify the strings (holes) of type j with 
the particles (antiparticles) of charge j {—j) of the Fermi gas. In other words, with each solution 
to the BAE for the finitized model corresponds a configuration of the Fermi gas, and vice versa. 
We remark here that we can in fact easily establish a bijection between Bethe Ansatz solutions 
and Coulomb gas configurations. So, for e.g., minimal configurations in the Fermi gas correspond 
to solutions without holes, and the leftward motion of a particle with charge j corresponds to 
shifting the Takahashi numbers by "inserting" numbers 1^ corresponding to the creation of 
holes. Using the map from the Fermi gas configurations to the configuration sums, we have in fact 
a bijection between admissible sequences of spins (<7i, . . . , cr m +i) and solutions to the Bethe Ansatz 
equations labelled by their Takahashi numbers. In this context we should mention that several 
authors |37], [Tj| [□], |2T| have indeed conjectured and/or computed the total number of solutions to 
the Bethe equations to be expression fl3.12j ), here counting the number of charge configurations in 
the Fermi gas. 

To further obtain the equations of TBA, we first let m — > oo in the Fermi gas partition function 
( |3.27|) , and then take the limit q — > 1~. This procedure, developed in refs. [[3^, |3{|, is equivalent to 
computing the T —>■ asymptotics of the entropy in TBA. Writing lim,™^^ 
have the asymptotics 



M) = Efcflfe? , we 



log a k 



1 



3tt 



S(T) 



T 



(5.6) 



T=0 



Computing a k by steepest descend gives us [15 



S(T) = -— 5>&) - L(Vj)\ + o(T), 

71 3=1 



(5.7) 



with L the Rogers dilogarithm function |40 
log(l-C) , logC 



Liz) 



2 Jo 



C 



i-C 



d(. 



(5A 



The numbers £j follow from the TBA equations 0, 
21og(l-^) = E^ log&), 



(5.9) 



k=l 



which arise here as the conditions for the saddle point. The rjj+i follow from the TBA equations 
with r replaced by r — 1, and rji = 1. 
Using the identity [(41 



r-3 



7T" 



1 - 



-5-1, 

r(r — 1)1 



(5.10) 



we finally get 



S(T) 




6 



t=o " v r[r - I) 



(5.11) 



in accordance with the TBA result of BR. From this last result the central charge of the ABF 

(5.12) 



models can be read off as 
6 



r(r — 1) 

To conclude this section, we wish to make a comment on the TBA calculations of BR as carried 
out in ref. ||. One of the shortcomings of their calculation is the failure to also derive results for 
the scaling dimensions of the ABF model. However, from the equivalence between the Fermi gas 
and the TBA approach, it is clear that the constraint equations ( |3.11| ) which follow from the string 
hypothesis (|5.2|) only describe the vacuum or groundstate sector of the model. As argued before, 
to compute other configuration sums than X m (l, 1,2), (which corresponds to computing scaling 
dimensions) additional boundary particles have to be introduced on top of the Fermi gas. These 
extra particles result in a modification of the constraint equations ( |3.11| ).p| We therefore believe that 



the string hypothesis as formulated by BR is incomplete and does not describe all solutions to the 
BAE. Indeed, the total number of eigenvalues of the row-to-row transfer matrix of the ABF model 
exceeds the number of configurations of the Fermi gas. The fractional structure of the additional 
boundary particles needed for the general configuration sums, seems to suggest the existence of 
half-strings, which have roots in the upper or lower part of the complex plane only.[] 



6 Summary and discussion 

In this paper we have presented a new method for computing the one-dimensional configuration 
sums of the ABF model in regime III. Our approach, based on the interpretation of the configuration 
sums as the grand-canonical partition function of a one-dimensional gas of charged fermions, results 
in polynomial expressions of fermionic type. This opposed to the bosonic polynomials obtained by 
employing the recurrence method of Andrews, Baxter and Forrester 0. Combining both bosonic 



and fermionic techniques leads to a proof of Melzer's polynomial identities ||10fl , different from 
the recursive proof by Berkovich |12[]. In taking the thermodynamic limit, this also proofs the 
Rogers-Ramanujan identities for unitary minimal Virasoro characters as conjectured by Kedem et 
al. JTT]. Since our fermionic method is mathematically equivalent to the thermodynamic Bethe 
Ansatz calculations for the ABF model of Bazhanov and Reshetikhin we have established a 
unification of the corner transfer matrix and of the TBA technique. Here it is to be noted that in 
our calculations the particle content of the model follows without making the usual string hypothesis 
characterizing TBA calculations. 

Although we have applied our method to the simplest possible configuration sum of the ABF 
model, the other sums can be computed be introducing additional boundary particles to the Fermi 
gas. The details of this procedure will be treated in part II of this paper What is more, our 



6 These modifications correspond to the "additional inhomogeneities" discussed in refs. 25 . 
7 In a recent preprint the generalization of TBA to compute scaling dimensions was also briefly mentioned. 
How it relates to the half-strings proposed here remains unclear at present. 



method is by no means restricted to the ABF model, and we believe that all restricted solid-on- 
solid models for which the TBA program has been carried out @, |5], || |45|, [46| admit a fermionic 
computation of their configuration sums. We thus hope that our approach leads to proofs of the 



many polynomial and Rogers-Ramanujan identities conjectured in the recent literature |L6L O, 24 



2£ 



We also note that our Fermi gas formulation of the ABF configuration sums can be reformulated 
to obtain a partition theoretic proof of Melzer's identities. Such a proof corresponds to establishing 
a bijection between partitions with prescribed hook differences |47] and charge configurations of the 
Fermi gas. Work on this bijection is currently in progress in collaboration with O. Foda [1~8"|. 



Finally, it would be extremely interesting to connect the approach to fermionic character rep- 
resentations of Bouwknegt et al. [49-51] using Yangian symmetries, and of Georgiev using vertex 
operators [0], to that of this paper. 
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Figure 2: The movement of a particle of charge k through a particle of charge j > k. According to 
rule R3 the configurations in figure c and c' are to be identified. The dashed lines are the baselines 
as described on page 8. We note that the procedure for drawing the baselines in fact rules out 
possibility c' thus avoiding overcounting of configurations. 
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Figure 3: A typical charge configuration in which some of the particles have formed charge com- 
plexes. 




Figure 4: The construction of the charge configuration of figure 3 out of a minimal configuration 
by carrying out the steps described under M1-M6. 




Figure 5: The drawing of baselines in a charge complex to determine its particle content. The 
small arrows indicate how to divide the j-th baseline into zeroth baselines of smaller complexes, 
(a) Drawing of the zeroth baseline, (b) Drawing of the first baseline, (c) Drawing of the second 
baseline, (d) Drawing of the third baseline. 



2j 2k 

Figure 6: Particles of charge j and k with labelling of begin-, end- and interior points. 




